old versions/plots.R

In SophiaDobko/assimReservoirs: This package assimilates reservoir volumes in Ceará, Northeast Brazil

#' Plot contributing basins
#'
#' This function plots the identified contributing basins of identBasinsGauges
#' @param list_BG output of identBasinsGauges
#' @export

plotBasins <- function(list_BG){
  plot(list_BG$catch$geometry, col = "white", main = "Reservoir with contributing basins")
  plot(list_BG$res$geometry, col = "cadetblue4", border = "cadetblue4", add = T)
}


#' Plot rain gauges
#'
#' This function plots the identified rain gauges of identBasinsGauges
#' @param list_BG output of identBasinsGauges
#' @param distGauges distance in km around the contributing basins to look for rain gauges as it was used in identBasinsGauges, defaults to 30
#' @export

plotGauges <- function(list_BG, distGauges = 30){
  plot(list_BG$catch_buffer, border = "green", main = paste("Basins with rain gauges within", distGauges, "km"))
  plot(list_BG$catch$geometry, add = T, border = "red", col = "white")
  plot(list_BG$gauges_catch$geometry, add = T)
}


#' Plot idw interpolation
#'
#' This function plots the interpolated precipitation in the contributing basins
#' @param list_BG output of identBasinsGauges
#' @param list_idw output of idwRain
#' @import manipulate
#' @import raster
#' @export

plotIDW <- function(list_BG, list_idw){

  dailyRain <- list_idw$dailyRain_table
  idwRaster <- list_idw$idwRaster

  par(mar = c(2,2,2,2), oma = c(1,1,1,1))
  manipulate(
    plot(idwRaster[[day]], main = names(idwRaster)[[day]]),
    day = slider(1,5))
}

#' Plot trmm of subbasins
#'
#' This function plots the averaged precipitation in all subbasins of the contributing basins
#' @param trmm_means output of trmmRain, geospatial dataframe with mean TRMM precipitation of each field of a given sf geometric object
#' @import sf
#' @export

plotTRMM <- function(trmm_means){

  par(oma = c(0,0,1,1))
  plot(trmm_means["trmm_mean"], border = "black", main = "Mean TRMM precipitation")
}

#' Plot strategic reservoirs
#'
#' This function plots the strategic reservoirs identified by resRouting
#' @param list_BG output of identBasinsGauges
#' @export

plotStratRes <- function(list_BG, list_routing){

  catch <- list_BG$catch
  res_main <- list_routing$res_main
  riv_catch <- list_routing$riv_catch
  reservoirs <- list_routing$reservoirs


  plot(catch$geometry, col = "white")
  plot(riv_catch$geometry, col = "cadetblue", add = T)
  plot(reservoirs$geometry, col = "gray55", border = "gray55", add = T)
  plot(reservoirs$geometry[reservoirs$id_jrc %in% res_main$id_jrc], col = "cadetblue4", border = "cadetblue4", add = T)
  legend("bottomright", inset = 0.01, fill = c("cadetblue4", "gray55"), legend = c("strategic reservoirs", "non-strategic reservoirs"), cex = 0.6)
}